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I. INTRODUCTION 

Theoretical models involving two orbitals per unit 
cell in one dimension have been employed to under- 
stand various strongly correlated systems e.g. heavy 
fcrmion compounds high-T c superconductors pHH 
spin-ladders p^Ell and Hubbard-laddersP^H^' These mod- 
els are often referred to as "ladder models" because 
one may visualize the two states as opposite sides of a 
rung. Recent experiments demonstrate the formation 
of self-assembled transition-metal nanowire s, add ing an- 
other possible realization of ladder models JHEU i n such 
systems, one-dimensional (ID) nanowires composed of 
transition-metal atoms are confined at step edges on sub- 
strate surface. In many physically relevant cases, the 
surface bandgap structure of the substrate material is 
such that the electronic states of the adatoms are decou- 
pled from the bulk substrate bands (at least to leading 
order and for low energies) and the <i-orbital bands de- 
rived from the transition-metal nanowire form a multi- 
component one-dimensional Fermi gas with variety of in- 
teractions. Classifying the kinds of behavior which may 
be observed in these systems is an important open ques- 
tion. In this paper, we will focus on the case where 
only two bands cross the Fermi energy. This situation 
is equivalent to two-leg ladder models, which have been 
previously studied. 

Past studies revealed that two-leg ladder sys- 
tems may ex hibit dynamical symmetry enlarge ment 
(DSE) 10 ' 2 ° ' 2 -%nd dualities among ground statesPEH 
DSE occurs when a RG flow leads to an effective low 
energy fixed point which exhibits a higher symmetry 
than that of original lattice Hamiltonian. In two-leg 
ladder models, the emergent symmetry is known to be 
0(6)xU(l) when the system is away from half filling, 
and 0(8) at half filling. The different ground states 
of these low-energy theories are related to each other 
by duality mappings, which are generalizations of the 
Kramers- Wannier duality seen in the two dimensional 
Ising model.^ 



We expect that many qualitative aspects of two-leg 
ladder models, including DSE and duality properties, 
hold also in our models of transition-metal nanowires. 
However, new features of transition-metal wires require 
additional investigation. The orbital degeneracy of the 
transition-metal d-levels permits a rich set of onsite 
Coulomb interactions. In particular, the Hund coupling 
favors locally high spin configurations, and we specu- 
late that the ground states may exhibit non-zero spin 
structure. Furthermore, the pair-hopping term, which 
has been neglected in some previous studies, is found 
to be important; it implies the system has a unique 
orbital symmetry. We generically expect that orbital 
symmetry breaking at the level of the one-body terms 
leads to velocity differences between different orbitals. 
A large velocity difference complicates the application 
of established methods such as refermionizatiorpS' to our 
model. Strong correlation may also lead to spin-polarized 
states or to ferromagnetism observed, for example, in Co 
nanowiresP^ Finally, the possible realization of an orbital 
selective Mott phase^in one-dimension is an open issue. 

In order to understand the effect of these fea- 
tures of transition-metal nanowires, we study the one- 
dimensional two-orbital Hubbard model using perturba- 
tive renormalization group (RG) and bosonization ap- 
proaches. We combine the ideas of DSE and duality to 
list the possible ground states of our models. We find the 
form of interaction relevant to transition-metal d-levels 
leads to a new group of 8 insulating phases when the 
two orbitals are completely degenerate. We then obtain 
ground state phase diagrams using RG and bosonization. 
In physically relevant parameter regimes, the stability 
of the ground states to velocity differences is also inves- 
tigated. To be complete, the fully spin-polarized case, 
where the model is reduced to a Hubbard model with or 
without magnetic field, is also briefly analyzed. 

The methods we employ in this paper are strictly ap- 
plicable when the system is truly one-dimensional, and 
the interaction is weak. Our results are complimentary 
to those we previously obtained for the same model using 
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mean- field approaches The present results allow more 
complete understanding of the weak-coupling regime, but 
the mean-field theory can treat the intermediate and 
strong coupling regimes. 

The organization of this paper is as follows. In Sec. [n] 
we explain the model and approximation we employed. 
In S ec. |III[ we derive bosonized forms of the model. Sec- 
tion [TVpTdevoted to explanation of various possible order 
parameters. In section [Vj we use the idea of dynamical 
symmetry enlargement and duality to understand the re- 
lations among ground states. Quantum phase transitions 
among these ground states are explained. In section |V*l| 
RG equations and obtained phase diagrams are shown. 
The fully-spin polarized case is briefly discussed in Sec. 
|VII| Finally, Sec. |VIII| is a conclusion and summary. 



II. MODEL 

We start from a multi-orbital Hubbard model repre- 
senting the transition-metal d-orbitals with local on-site 
Coulomb interactions 



H 



EE 

m,s 



h.c. 



(1) 



Here c^ s is the annihilation (creation) operator for a d- 
electron in orbital m with spin s at site i. t™ m is the 
hopping between from orbital m on site i to orbital ml 
on site j. The interaction terms iJ; n t will be shown be- 
low. We believe this model encapsulates the physics of 
transition-metal nanowires. Through the paper, we set 
the lattice constant to 1. A symmetry breaking field oc- 
curs due to the one-dimensional geometry, and the pres- 
ence of the substrate may further lift the degeneracy of 
orbitals as well as providing an arbitrary ionization level. 
Thus, in a general case, one may have many d-derived 
bands with an arbitrary Fermi energy. 

For the sake of simplicity, we will consider here only 
the cases where two orbitals, A and B, are present at 
the Fermi level. The rotational symmetry in H lrA as we 
will see always allows us to diagonalize the hopping ma- 
trix, so we will ignore off-diagonal terms in the hopping. 
The band structure is then characterized by four Fermi 
points: two Fermi momenta, kji and fcs, and two chiral- 
ities r = R, L, representing electrons around the positive 
(i?) and negative (L) Fermi momentum. The total par- 
ticle number is n = 2(fc J 4 + ks)/^- In principle there 
are 5 possible cases which are summarized in Table [TJ In 
cases (a) and (b), the two Fermi momenta are equal, and 
the filling is commensurate and incommensurate respec- 
tively. In cases (c) and (d), the two Fermi momenta are 
different, while (c) is at half filling, and (d) is away from 
half filling. Finally in case (e) , one band is commensurate 
and the other is not, allowing an orbital selective Mott 
state. 

For the two-orbital system the interaction terms have 



the following form: 



i,m 

+ U'^n iAs n lB s 

i,s 

+ (V - J) n iAs niBs 

i,s 

i,s, 



(2) 



where and 



is the electron density. 



s = — s and m means B(A) if m = A(B). Since the 
substrate will screen the long-ranged Coulomb interac- 
tion, we include only on-site Coulomb interactions. U 
and U' indicates on-site Coulomb repulsion between two 
electrons in the same band or different bands, and J rep- 
resents the Hund coupling favoring high spin state. J' is 
the so called pair-hopping term. For a transition-metal 
ion in free space all of these parameters are positive. 
Now, we assume the following relationship among them, 
which is usually preserved among d-orbitals, 



J' = J 
U' = U- 2 J. 



(3) 
(4) 



To make the symmetry of the interaction terms ex- 
plicite, we introduce the following charge, spin, and or- 
bital (pseudo-spin) operators: 



ms 

Si = — Pirns' 7 ss'Cims' , 

mss' 

Ms 



^m^mm' C-irn' s j 



(5) 
(6) 
(7) 



where er and r are Pauli matrices. Then, the interactions 
in terms of U and J are given by, 



U 



E + JS 



- 1 3JTf 



2J(Tff 



U + 5J 



(8) 



Since U' > 0, the physically relevant range of J is limited 
to < J < U/2. 

Now we discuss the symmetry of our Hamiltonian at 
the bare level. When J = 0, the interaction term pos- 
sesses U(l) c xSU(4) S!0 x Z 2 ; the indices "c", "s", and 
"o" denote the charge, spin and orbital parts, and the 
Zi symmetry refers to the interchange of orbitals. When 
J ^ 0, the symmetry of the spin-orbital part is broken 
to SU(2) S xU(l) D . We'd like to emphasize that in our 
convention, the U(l) axis is the orbital y-axis, wheres 
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Case 


Fermi momentum 




symmetries 


(a) 


kA = k B 


= 2 


U(l) c xSU(2) s xU(l) x2 2 


(b) 


k A = k B 


/ 2 


U(l) fl xU(l)LxSU(2) s xU(l) x2 2 


(c) 


k A / k B 


= 2 


U(l) c xSU(2) fl xU(T) o 


(d) 


k A / k B 




U(l)«xU(l)z,xSU(2).xU(l) 


(e) 


kA / k B = 7r/2 


/2 


U(l) c xSU(2) s xU(f) o 



TABLE I. Possible band structures in the two-orbital Hubbard model, and its symmetries when va = v B . U(l) r represents a 
gauge transformation of particles with chirality r. 



B A 



A B 









FIG. 1. Various scattering processes and "g-ology" 



it is orbital z-axis in the convention of Refs. [33] and 
1291 Only if va = vb & n d kA — ks, as in cases (a) and 
(b) , does the total Hamiltonian have the same symmetry 
as the interaction. Otherwise, including cases (c)-(e), 
the total symmetry is reduced to U(l) c xSU(2) s at the 
bare level due to the lower symmetry of the kinetic term. 
However, it is known that low-energy theories in weak- 
coupling still have an effective U(1) Q symmetr ^ 10 * 21 * ^ at 
least ii va = vb- We will see this in detail later. 



III. BOSONIZATION 

In ID systems, it is known that bosonization enables us 
to describe the low energy physics in a simple manner. In 
this section we present an Abelian-bosonization analysis, 
which is useful for the cases when difference in velocities 
is negligible. Appendix [B] outlines the more complicated 
formalism needed for unequal velocities. 

To classify the various scattering processes we follow 
the notation of Ref. "1-4" corresponds to conven- 

tional "g-ology" indices for left and right moving fields, 
and "<z-d" are similar but label orbital indices. Some 
examples are given in Fig. [TJ The SU(2) S symmetry 
constrains coupling constants as 



m _ m _ n 
Sid 92d — u 

92b - 92b ~9la = 

9ic-92c+9\\ c = 0. 



(9) 



For the rest of the paper, we will omit _L when it is not 
confusing. 



When two Fermi points coincide, we have additional 
processes gii a , gib, and g2a, which are connected by the 
SU(2) S symmetry as, 



.911 a - 9lb + 92a = 0. 



(10) 



f^A = fc-B means that we have the U(1) D symmetry, and 
this implies 



~92d 
9 Id 



9l(2)d - 9l(2)d — 
H 92b + 92c + 92a = 

H gib + gic + gia = 

92b +9\\c - 9\\a= 0. 



(11) 



As Umklapp processes, we have gnu and g^ (i — a- 
d). The g™ d process represents the intraband Umklapp 
process, so it only exists for k m = tt/2. Other Umk- 
lapp processes are possible whenever n — 2. The SU(2) S 
symmetry gives 

.93a +9\\b -.936 = 0, (12) 
and the U(1) D symmetry leaves 

93d = 93a + 93b + 93c- (13) 

Using standard bosonizationp^H ^ ne Hamiltonian 
density Ho of the free-boson part becomes, 

1 



- y 



2ir ^ V » v 

fJi — C.S 
V=Q,7T 



K„ 



(14) 



where <fi an d are connected to density and current: 
oc n, and V# oc j. The renormalized Luttinger pa- 
rameters and velocities are given by 



K. 



cO(tt) 



V c 0(tt) 
v s0{-k) 





(yid 


± y 26 )/2 


1 + 


(yid 


±V2b)/2 


1+ 


(yid 


±lft«)/2 




{yid 


±Vla)/2 



= 1 



VcOi-n) 



J/s0(7r) 



(15) 



Vi - {yid ± V2b) 2 /4 
V 1 - (yidTyia) 2 / 4 



with y, = gi/4nv and y 2b = y^ b + y\ b . "c" and "s" rep- 
resent the charge and spin modes, and "0" and 'V are 
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used for bonding and antibonding combination. We in- 
troduced y^v for each Luttinger parameter for later use. 
A detailed derivation of this Hamiltonian is given in Ap- 
pendix [AJ We ignore the velocity difference induced by 
the <?4 process, since its effect is to shift the phase bound- 
aries slightly. We will provide a separate treatment for 
systems in which difference of initial velocities is quite 
large. 

The SU(2) S symmetry fixes K s0 ^ to be 

Kso^) = 1 + (Vid ± yia)/2 = 1 - y s0 (Tr) (16) 

along the RG flow. Similarly the U(l) c symmetry (when 
it exists) constrains K C7r to be 



Ken = 1 + (-yu + V2b)/2 = l-y c 



(17) 



The interaction part of the Hamiltonian is rather com- 
plicated. The interaction terms common to all cases are 



Him = gid J cos (20 sO ) cos (2(j> sw ) 
+ 9ia J cos (2</> s0 ) cos (20 S7r ) 

- 9xc J cos(20 sO )cos(26» CTr ) (18) 

- 52c J cos (20 S7r ) cos (20 C7r ) 
+ 9\\c / cos (26» S7r ) cos (20 C7r ) . 



Here J stands for ^a) 2 I an< ^ a ^ s ^ ne cut-off which 
is of the order of the lattice constant. The last term 
does not exist in the original Hamiltonian, but will be 
generated after renormalization. 

When k,A — ks, in case (a) and (b), additional pro- 
cesses are allowed, 



H int = 9\\a J COS {2<j) sw ) COS (2</> c7r ) 

+ J cos(2(/> sO )cos(20 CTr ) (19) 
+ g2a J cos (26» S7r )cos(20 c;r ). 
Finally when the filling is commensurate, we have 
H- nt = g 3a J cos {2<j) cQ ) cos {26 S7I ) 
~93b J cos (2</> c0 ) cos (2^ S7r ) 
-53c / cos (2^ c0 ) cos (20 C7r ) (20) 
~.93d J cos (20 cO ) cos (20 OT ) 
-g\\ b / cos (20 cO ) cos (2</> s0 ) ■ 



The 33d process exists only when both bands are com- 
mensurate, case (a). Again, we ignored the type inter- 
actions whose scaling dimension is always larger than 2, 
and this is consistent with the equal velocity approxima- 
tion. The initial values of these coupling constants are: 

gid = 93d = 4U, gia = gic = 52c = 52a = 53a = 53c = 4 J, 

526 = 9\\a = 9\\b = 4(C/ - 3J), and g u = g^ b = g 3b = 
4(U — 2J). They will take different values after renor- 
malization. 

Finally we will discuss the symmetry of the linearized 
model (Table [j]!. First, U(l) c and SU(2) S (around z-axis) 
are displayed in the invariance of the Hamiltonian under 
the translation of 9 c q and 9 s q. In fermionic language, 
each corresponds to the following gauge transformation: 



->■ e 



(21) 



The indices r, to, and s represent chirality, orbital, and 
spin, and a expresses the constant phase shift. The con- 
served Noether currents, corresponding to the U(l) c and 
SU(2) S symmetries are, 



N rms oc J dxV<j) c0 

rms 

^ sN rms oc J dxV<j) S Q 



(22) 



where N rms is the particle number at branch specified by 
r, m, and s. Since V(/> is momentum conjugate of 6, the 
operator exp (J dxVcj)) gives a constant shift of 8. The 
SU(2) S rotation around x- and y-axis are not manifest in 
Abelian bosonization. 

Away from half filling, there is also a continuous chiral 
symmetry under the transformation, c rms — > e lra c rms . 
Thus, the Hamiltonian is invariant under arbitrary trans- 
lation of c o, with conserved total currents, J a + Jb oc 
/ dxV9 c0 where J m = J2 S N Rms ~ N Lms . At half filling, 
this symmetry is broken to a discrete symmetry, and true 
long-range order can be realized.^ When a system has 
both the chiral symmetry and the U(l) c symmetry, this 
implies that left, and right moving parts have separate 
conservation laws corresponding to the U(l)ijxU(l)i 
symmetry. Similarly, other gauge transformations such 
as 



irsa 

^rms 7 c ^rms 



. irmsa 
^rras 7 l - L rms 



(23) 



leave the Hamiltonian invariant for discrete values of a 
mr/2; each corresponds to the discrete shift of 9 cn , 6 S 
4>s0, and S7r . 



The effective 



ks and va = Vb 



U(l)„ sym metry appearing when Ua 7^ 
d EHElHlSl corresponds to the invariance 
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under the translation of </> C7r or the gauge transformation, 

irma /94 s ) 
^rms 7 c ^rms- \ ) 

The conserved "charge" corresponding to this symme- 
try is the difference of two orbital currents: J a — Jb oc 
/ dx\78 C7r . 

When kA = fcs, there is an explicit orbital rotational 
symmetry about y-axis. This transformation mixes 
fermions in different orbitals, so its generator cannot be 
expressed as a local operator in Abelian bosonization. 
This fact leads to a new combination of possible ground 
states as will be shown. 



IV. ORDER PARAMETERS 

A. Order parameters away from half filling 

As order parameters, we take fermion bilinears char- 
acterized by the chirality, spin and orbital indices; thus 
in the model considered here there are particle-hole bi- 
linears, 

(Aph) rr < = c[. ms C r i m i s ' (25) 

and particle-particle bilinears, 

( A Pp)rr' ;I ™ m = msc l ms cl,^rjr, (26) 

where elms is the annihilation (creation) operator of elec- 
tron with chirality r, orbital m, and spin s. Since combi- 
nations of r = r' are irrelevant in a RG sense, we will fix 
(r, r') — (R, L) in the remainder of this paper. We will 
use the following convenient basis to represent them, 

0%= E T* m X S '(Ap h r W +h.c. (27) 

mm' ss' 

Og = X) r^ mt ai s ,(\ P r' imm ' +h.c. (28) 

mm' ss' 

where i,j — (0,1,2,3) and r and a are Pauli matrices 
with T° b = cr° & = 6 a b- These transform as rank 2 ten- 
sors under SO(4) ~ SU(2) s xSU(2) transformations; the 
SU(2) S rotations connect a 1 ' 2 ' 3 , and the U(1) D rotation 
(if it exists) connects r 1 and t 3 . Thus, we take the quan- 
tization axis along z-direction for spins. Additionally, 
as we will show below, all the high spin states (j = 3) 
such as SDW states and triplet superconductivities are 
excluded from the possible ground states in our models, 
so we will not consider them. 

We label our order parameters by the transferred mo- 
mentum and the sign at each Fermi point. Here we have 
4 Fermi points each degenerate about spins, so, in princi- 
ple, there are 4 possible cases (Table |llj) . We use s-wave 
when all 4 points have the same signs. p x and p y are odd 
under the inversion R O L, and A O B respectively. 
e?-wave is odd under both inversions. 



Applying this classification, we find that i = 0, 1 are 
both s-wave for the particle-hole channel, while the for- 
mer is intraband type and the latter is interband type. 
We put " ' " for an interband order to distinguish these 
two. i = 2 is found to be interband py-w&ve, and i = 3 
is intraband p y -wsve. For the particle-particle channel, 
we found d'-, p y -, s-, and s'-wave orders for i = 0, . . . , 3 
accordingly. Since p^-wave does not appear, we will use 
p for p^-wave orders in the rest of the paper. We note 
that the d-wave superconductivity, which often appears 
in two-leg ladder models^ is the pSS state in our nota- 
tion. The result is summarized in Table [TTT1 

The order parameters in bosonized forms are, 

Oph cos ( k ox - 4>co) sin (k w x - </> C7r ) sin (0 sO ) sin (0 S7r ) 

OpS oc cos (k x - (f) c0 ) sin (6 cn ) cos (<f> s0 ) cos (0 S7r ) 

OpS oc cos (k x - (t> c0 ) cos (0 C7r ) cos ((j> s0 ) cos (0 s7r ) 

Oph oc cos (k x - <j> c0 ) sin (k^x - $6 C7r ) cos (^ s0 ) cos (0 S7r ) 

(29) 

0°p oc e~ l9 "° sin (k % x - <p C7I ) cos (0 sO ) cos (0 S7r ) 

OH « e ~ l8c ° sin ( on) sin (<M sin (4> sn ) 

2 p ° p K e ~ l9c " cos sin (^o) sin (0 S7r ) ' '"' ' ' 

oH K e ' i0cO sin ( k * x ~ ^ sin (^°) sin ' 

with fco(,r) = kA ± kg- If coupling constants grow to 
the order of t after renormalization, the corresponding 
bosonic fields are pinned to the values which minimize 
the resultant potential. For incommensurate fillings, 
the total charge mode is massless, and the interaction 
terms pin the other modes to definite values (in mod of 
7r/2). The pinned values determine the order parameter 
which gives a finite value of correlation. It is easy to 
find such pinned values from the above expressions e.g. 
(0 c ^,(j) s o,(t) S7T ) = (71-/2,0,0) for the s'CDW state. 

B. Order parameters at half filling 

When a band is commensurate, the number of possible 
ground states increases since the total charge mode also 
comes into play. After renormalization, all modes become 
massive unless the system remains in a Luttinger Liquid 
phase. The phases are classified by the fixed values of 
pinned fields. 



"Angular momentum" 


(A,R) 


(A,L) (B,R) (B, L) 


s 


+ 


+ + + 


Px 


+ 


- + - 


Py 


+ 


+ 


d 


+ 


+ 



TABLE II. Angular momentum and sign distributions 



G 





Types of particle-hole order 


Types of particle-particle order 


(0,0) 


Charge density wave (CDW) 


d'-wave singlet SC(d'SS) 


(1,0) 


s'-wave charge density wave (s'CDW) 


Py-wave singlet SC (pSS) 


(2,0) 


p^-wave charge density wave (p'CDW) 


s-wave singlet SC (sSS) 


(3,0) 


Py-wave charge density wave (pCDW) 


s'-wave singlet SC (s'SS) 



TABLE III. Classification of order parameters. " ' " indicates that the order is inter-orbital type. We use p to express p^-wave 
orders, since p^-wave orders do not appear. 



The first set of new order parameters we consider is 
"bond" orders, which is basically the density-wave slid 
from on-site to "on-bond" . Away from half filling, site 
and bond density waves are degenerate and the total 
charge mode is massless. However, at half filling, these 
two states decouple, and their expectation values of (</> c o) 
arc different by tt/2. We call these on-bond states BDW 
(bond density wave). Each phase is doubly degenerate 
about the translation of (4> c o) by 7r. 



and B at site i, we find 

|«SS) = Y[ [CiA-fCiAl + CiBfCiBl] |0) 
i 

\pSS) = Y[ [CiAfCiAi. - CiB-fCiBl] |0) . 
i 



A second class of insulating phases evolves into super- 
conducting phases upon doping. These insulating phases 
share similar properties with the superconducting phases 
in the sense that the pinned fields and values are the 
same except for the total charge mode. While we will 
explain each state in detail later, we show the results 
first. At half filling, the s'SS state becomes the Haldane 
charge (HC^state or the rung-singlet (RS) stateP^Sim- 
ilarly, the d'SS state turns into the Haldane-orbital (HO) 
stateP or the rung-triplet (RT) stateP The S-Mott and 
S'-Mott phases are developed from the sSS state, and the 
D-Mott and D'-Mott phases are from the pSS state. Two 
insulating phases developed from the same superconduct- 
ing state are connected by the shift of 4> c q by tt/2. 

Another way of looking at these insulating phases is 
to regard them as Ising disorder phases of the CDW and 
BDW phasesP The disordered states can be obtained 
by applying the mapping, cj) S7T f-> 6* S7r , to density- wave 
states. While the order parameters of density- wave states 
can be expressed by the local fermions, these disordered 
phases need non-local string operators. The Ising dual 
phases of the s'CDW(BDW) state is the S(')-Mott state, 
and the p'CDW(BDW) state is dual to the D(')-Mott 
state. The RS, RT, HC, and HO states are Ising dual 
to the BDW, pBDW, CDW, and pCDW states respec- 
tively. These disordered states are non-degenerate un- 
like density-wave ground states. The correspondence be- 
tween insulating/metallic phases and ordered/disordered 
phases is summarized in Table [TV] The pinned values of 
16 insulating phases can be found in Ref. [231 



These two-particle states are the same as the ones given 
in Ref. [16] as the S-Mott and D-Mott states in the bond- 
ing/antibonding basis. Therefore, we identify the cor- 
responding insulating states as the S-Mott and D-Mott 
states. The S'-Mott and D'-Mott states are different from 
these only by 64> c q — 7r/2, and this indicates that the 
order is located on bond. To be precise, while the au- 
thors of Ref. HI] claimed that the D-Mott and S'-Mott 
states are connected to the RS and RT states in the 
strong coupling limit, this is not the case in our mod- 
els as we will see later. The discrepancy comes from 
the fact that the strong transverse hopping requires the 
use of bonding/antibonding basis in weak coupling anal- 
ysis, although the strong coupling analysis starts from 
two independent ladders. However, in the weak coupling 
regime, the properties of Mott states are the same. 

Next, the s'SS and d'SS states can be represented by 



Ordered 


Disrdered 


Insulating 


Metallic 


Insulating 


Metallic 


CDW 
BDW(SP) 


CDW 


HC 
RS 


s'SS 


s'CDW 
s'BDW 


s'CDW 


S-Mott 
S'-Mott 


sSS 


p'CDW(SF) 
p'BDW (FDW) 


p'CDW 


D-Mott 
D'-Mott 


pSS 


pCDW 
pBDW(SP^) 


pCDW 


HO 
RT 


d'SS 



To see the above properties of each phase more care- 
fully, we consider two-particle states corresponding to 
each superconducting state. We first look at the sSS 
and pSS states. Considering the occupation of orbital A 



TABLE IV. Correspondence between insulating/metallic 
phases and ordered/disordered phases. The commonly used 
names are given in parenthesis. SP: spin-Peierls, SF: stag- 
gered flux, FDW: /-wave density wave. 
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the following two particle states, 

|s'SS) = Y\_ [ c iA^CiBl - CiAi,C iB f] |0) 
i 

\d'SS) = Y]_ [CiA\C i+ x t Bl - CiAlCi +lf Bt - (A -H- B)} |0) , 

(32) 

where we ignored the total charge mode, and assumed 
hj\ = k B = 7r/2 for the d'SS state. Obviously, the s'SS 
state forms a singlet Cooper pair on a rung, while the 
d'SS state forms such singlet pair in staggered way. These 
pairings are the reminiscent of resonance-valence-bond 
(RVB) in the RS and RT st ates which appear in spin- 
\ Heisenberg two-leg laddersP^ES The RS state appears 
in chains with antiferromagnetic coupling along rung (or 
ferromagnetic coupling over plaquette diagonals), and it 
is a RVB state whose stable configurations are the singlet 
along rung or ladder as in Fig. |2(a) We can see the 



clear connection between the RS state and the s'SS states 
in Eq. ( 32 ) . This state is characterized by a non-zero 



expectation value of a string operator 

((Sl,i + S B>i+1 ) e-EfcK^+^O (<?* + S BJ+1 )), 

(33) 

and exhibits a "hidden antiferromagnetic order" ; as we 
group Sa,i and S B ,i+i, the total spin align antiferromag- 
netically along ladder except spin-0 site. 

For chains coupled ferromagnetically along rung (or 
antiferromagnetically over plaquette diagonals), the 
above order disappears, and the Valence-Bond -Solid 
(VBS) configuration^ becomes stable as in Fig. 2(b)| 



As can be seen, singlet pairs are formed in a staggered 
manner as in the d'SS state, and this results in triplet 
pair along a rung. Thus, we call it rung-triplet (RT). The 
following order parameter takes a non- vanishing value, 



((s z A 



S B>i ) e«S:l(5, t +«.) (S A<j + S BJ )), (34) 



which represents a hidden order about the spin-triplet 
on a rung [Fig. 2(b)| . This state smoothly continues to 
the Haldane gapped state as interchain coupling becomes 
large. 

The string operators for these disordered states are 
non-local, so complications arise in the transcription to 
bosonic variabl es. A generally accepted form for Eqs. 
(pi and pi} ifPEH 



(sin(0 sO (x)) sin(0 sO (2/))). 



(35) 



This correlation function takes a non-zero value also for 
the s'SS and d'SS phases. Considering these properties, 
we conclude that the RS and RT states are the corre- 
sponding insulating states of the s'SS and d'SS states. 
In particular, we found that the d'SS and RT states are 
the disordered versions of high spin states. As we men- 
tioned, the D-Mott and S'-Mott states do not evolve to 
the RS and RT states. 




-1 10-1 1-1 1 

(a) rung-singlet (RS) 

i f x |3 OK) f x i3 

cM €3 i 0-4 i 

-10 1-11-110 0-1 

(b) rung-triplet (RT) 

FIG. 2. Rung singlet and rung triplet states (after Ref. 133 j> . 
Closed circles represent singlet bonds, and dashed lines indi- 
cate hidden antiferromagnetic order. 



The bonding counterparts of the RS and RT phases 
are the Haldane-charge (HC), and Haldane-orbital (HO) 
phases proposed by Nonne et cdP3 as the Haldane gapped 
states of pseudo spin-1 antiferromagnetic Heisenberg 
chain; this realizes when the charge or orbital symme- 
try is promoted from U(l) to SU(2). The form of string 
operators for these states is similar to Eqs. (33 1 and (34) 



if we replace SU(2) spin operators by the charge and or- 
bital SU(2) operators, 



j; 



-(n A ,i + n B ,i), J- = 4.aA.bi ~ c 1al c L 



^i,A\^i,Bl ~ ^i,A±^i,B^ (36) 

T i = n A,i - n B ,i, T} = ct )j4t Cj,fl t + c\ Ai c t .Bi- (37) 
The bosonized form of string operators becomes 



(sin 
(sin 



.(0co(a;))sin(0 c o(y))) 
(x)) sin(0 

I'll (y))}, 



(38) 
(39) 



where V0 c o ~ J z (x) and V<^ C 7r ~ T z (x). These expres- 
sions agree with the fact that in the insulating phase 
the total charge mode should be massive, and that <f> C7T 
is fixed also in the s'SS and d'SS states. Further study 
of the connection between phases appearing in strong- 
coupling limit, and those in the weak coupling limit 
would be of interest, but will not be pursued here, be- 
cause our focus is on the weak coupling limit. 



DUALITY, AND DYNAMICAL SYMMETRY 
ENLARGEMENT 



A. Possible ground states 

Two of the most p romin ent features of one-dimensional 
systems are duality 22 23 and dynamical symmetry en- 
largement (DSE)^^The idea of duality is based on the 
observation that the low energy theory is invariant un- 
der some discrete operations apart from the continuous 
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n c 


n a 




0/ 




n s 




(a) 


BDW 


CDW 


p'CDW 


p'BDW 


S-Mott 


S'-Mott 


RT 


HO 


(b) 


CDW 


CDW 


p'CDW 


p'CDW 


sSS 


sSS 


d'SS 


d'SS 


(c) 


S-Mott 


S'-Mott 


D-Mott 


D'-Mott 


s'CDW 


s'BDW 


p'BDW 


p'CDW 


(d) 


sSS 


sSS 


pSS 


pSS 


s'CDW 


s'CDW 


p'CDW 


p'CDW 



TABLE V. Possible ground states for case (a)-(d) in Table |T| at equal velocities. The BDW, CDW, S-Mott, and sSS phases 
in the second column are expressed by the fundamental SO (8) or SO (6) Gross-Neveu model for case (a)-(d) respectively. The 
other states are mapped from these fundamental states by the duality transformation fl in the top row. Q„ is an operation of 
£1 — > — £™ for all the a's in a symmetry sector v. For example, in case (a), the BDW state is mapped to the p'CDW state by 
Q - For doped cases, (b) and (d), the charge mode is separated from the rest, and the ground states are invariant under Q c . 



symmetries listed in Table [T] These discrete symmetries 
enable us to relate one ground state to another, and to 
understand quantum phase transitions among them. 

DSE means that the effective theory describing the low 
energy fixed point exhibits a higher symmetry than that 
of the original lattice Hamiltonian. This phenomenon 
was noted by Lin et a£p^, who found that the low-energy 
theory of half-filled two-leg Hubbard ladder is the S0(8) 
Gross-Neveu (GN) model. Since their work , DSE has 
been seen in other multiband systems? 36 * 37 * Combining 
these two ideas, we will identify the possible ground 
states for our models, taking DSE for granted. 

Now, as a preparation, in order to exhibit the symme- 
tries of the Hamiltonian, we refermionize the model using 
8 Major ana fermions as explained in Refs. H31 and 
We decompose each mode into two Majorana fermions as 



i 



(40) 



At half filling, the obtained expression for the U(l) sym- 
metric case, (a), is, 



a=l 

9l 2 

+ '—n s + g 2 K s K, + g3K s m (41) 
, . 95 2 . 96 2 

+ + y k o + y K c 

+ g 7 n s n c + g%K K c + g 9 K!K c 



where k s = EL^Ul *o = TLitteh «/ = 
and k c = Ea=7^£' The indices V, "o", and "c" 
refer to the SU(2) S , U(l) , Z 2 , and U(l) c symmetries 
respectively. 

Away from half filling, case (b), the charge mode is 
decoupled from the other modes, and we don't need to 



consider gi = e^g. Thus, we have 



a=l 



9l 2 

1 . 95 2 

+ g i K K I + yK . 



93H s Kl 



(42) 



As we mentioned before, even without an explicit U(l) 
orbital symmetry for the lattice Hamiltonian, the low 

energy theory has the effective U(1) Q symmetry ! 10 * 21 ! 23 ! 
Therefore the structure of the refermionized forms are 
the same as above. Thus, cases (c) and (d) now have the 
same form as Eqs. (41) and (42) with different values of 
5 'sP 

For the refermionized forms, duality mappings are de- 
fined as, — > — while keeping right-moving parts 



untouched. It is easy to see that Eqs. (41) and (42) 



are invariant under such transformations if we change 
the signs of some coupling constants as well. To retain 
the form of the Hamiltonian, only mappings that trans- 
form all the Majorana fields in the same symmetry sector 
are permitted. For example, for the SU(2) spin part, we 
should map the three left Majorana fermions, ££ =1 ~ 3 , at 
the same time. For notational convenience, we define £l v 
as an operation of — > — for all the a's in a symmetry 
sector v. With this at hand, it is obvious that allowed 
mappings for half filling cases are 



0(8) 



= {^c, ^o, SI/, ^si ^c,o, ^c,J> ^o,l} ■ 



(43) 



The number of independent mappings is 3, and other 
mappings just follow from them e.g. fl j = f2 Sl/. Away 
from half filling, the charge mode is separated, so only 3 
of the above mappings are left, 



^0(6) = {^o, 0j, Cl s }, 



(44) 



and two of them are independent. An immediate conse- 
quence of these dualities and DSE is that although we 
showed 16 insulating phases for half filling systems, and 
8 metallic phases for incommensurate filling, only a part 
of them are realized. 

Now, we will show such possible ground states for each 
model. We start from the "fundamental" SO(8) GN 
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model, 



(45) 



which appears at low-energy when all the c/'s in Eq. (41 ) 



converge to the same value as a result of DSE. For case 
(a), this model represents the BDW phase, and other 
possible phases are found by applying Qo(8) to the BDW 
state (see Table |v| . We denote them as T y , and they are 

T y : BDW, CDW, p'CDW, p'BDW, 

S-Mott, S'-Mott, RT, HO. (46) 

The case (b) follows from the relation between insulating 
states and metallic states (Table 
to the CDW phase, which is "fundi 



IV|), or applying fi 0(6) 
: amental" . The original 
lattice model we are considering here is invariant under 
the orbital U(l) rotation about y-axis. As we mentioned 
in Sec. |III| the generator of this symmetry cannot be 
expressed by a single local bosonic field within Abelian- 
bosonization scheme. 

This combination, T y , is different from the ones which 
have been studied extensively; previously studied phases 
are 

T z : BDW, CDW, pBDW, pCDW, RS, HC, RT, HO 

(47) 

and 

f\J : S-Mott, S'-Mott, D-Mott, D'-Mott, 

s'CDW, s'BDW, p'BDW, p'CDW. (48) 

The former, T Z1 appears in models with weak trans verse 
hopping, and with the U(1) Q symmetry about z-axisP^HH 
The latter, T z , appears when the model has strong trans- 
verse hopping, and the low energy theory possesses the 

U(1) Q symmetry) 7 ! 10 ! 12 ! 16 ! 17 ! our cases (c) and (d) belong 
to this category (see Table |V| . 

The connection between T y and T z is obvious. Since 
the generator of the orbital symmetry for each case is y- 
or z-componcnt of Eq. ( 37 1 , they are simply mapped to 



each other by a rotation around x-axis: 



-rAs 



-rBs , 




(49) 



This transformation does not affect the charge and spin 
generators. For instance, the S'-Mott state in T y goes to 
the HC state in T z by R x . The correspondence among 
other states is given in Fig. [ij On the other hand, T z and 
F z transfo rm ea ch other by so-called strong-weak tunnel- 
ing duality^EQ 







CLmt -> C' LmV C Lml -> -C' Lmt . (50) 

The relationships among metallic ground states are easily 
obtained by Table \N\ 

Therefore, we found that in addition to underlying 
band structure, the form of interaction also affects the 
possible combinations of ground states. We summarized 
these results in Table [V] and Fig. [3] 




BDW, CDW, p'CDW, p'BDW, 
S-Mott, S'-Mott, RT, HO 




R x and Q I 



BDW, CDW, pBDW, pCDW, 
RS, HC, RT, HO 




S-Mott, S'-Mott, D-Mott, D'-Mott, 
s'CDW, s'BDW, p'BDW, p'CDW 



FIG. 3. Relationships among three groups of insulating 
ground states. R x indicates the rotation about x-axis in or- 
bital space in Eq. ( 49 1 , a nd Q± is strong- weak tunneling 
duality mapping in Eq. (50 1. For example, the S'-Mott state 
in r y maps to the HC state in F z by R x . 



B. Quantum phase transition 

The quantum phase transitions among gapped ground 
states could be either first order or second order. For the 
transitions among states connected by duality, the modes 
which are not involved in the mapping become massive at 
higher energy, and the effective low energy theory near 
the transition contains only Majorana fields flipped by 
the mappingPsl 

For a single Majorana field, it becomes the critical Ising 
model, 

H = (ZrMb. ~ - imtHti. (51) 

Over the transition, the mass changes its sign, and this 
represents a second-order phase transition. With more 
than one field, the low-energy effective theory becomes 
the massive O(N) GN model, 



(52) 

The fate of further renormalization to lower energy de- 
termines whether the phase transition is first-order or 
second-order depending on the final fixed point for the 
critical fields! 10 1 16 ' 21 1 39 ' ^ 1 The transition line is defined as 
the point where the m in Eq. [52] goes to zero, and the 
critical fields are expressed by a massless GN model in 
the vicinity of transition. For N = 2, it is known that 
the system can be mapped to a Gaussian model, so it is a 
second order transition. For N > 3, however, if the cou- 
pling constant in the GN model is positive (g > 0), the 
renormalization flow departs to a strong coupling fixed 
point (asymptotic free), since RG equation is given by 



9 oc 9 



(53) 
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At this fixed point, the mass is generated dynamically, 
and the system is off-critical. We can see this either 
by mean-field treatment of the interaction (reducing the 
quartic part to quadratic with the order parameter, 
(£h£l))j or by stationary phase approximation, which 
becomes exact when N — > oo. At this massive fixed 
point, there are two degenerate minima about two signs 
of mass, and they correspond to two phases connected 
by this first-order transition. On the other hand, when 
g < 0, further renormalization reduces g to 0, and the 
system reaches a massless fixed point; this represents a 
second-order transition. 



When the transition is second order, the critical theory 
is described by a conformal field theory (CFT) due to its 
dimensionality, (1+1). Each CFT is characterized by its 
central charge, c, which roughly expresses the number of 
critical fields, c = 1/2 is the Ising critical theory, c = 1 
is the U(l) Gaussian theory, and c = 3/2 is the SU(2) 2 
Wess-Zumino-Novikov-Witten theory. With the duality 
mappings, it is easy to read off the central charge of each 
CFT. Since each Majorana fermion carries c = 1/2, the 
number of fields flipped by a mapping directly tells us 
the central charge. We will identify the phase transitions 
appearing in our phase diagrams more precisely in the 



next section. 



VI. RG EQUATIONS AND PHASE DIAGRAMS 

In the RG equations obtained from Abelian bosoniza- 
tion, we use normalized coupling constants defined as 



y, 



9i 

4ttv 



(54) 



RG equations ar e deriv ed using the operator product 
expansion (OPE^^HUiSI anc j integrating out higher fre- 
quency modes. The RG equations are complicated, so we 
will not show them explicitly. In all of the cases which 
we have examined, the RG equations may be expressed 



dm 
dl 



dV 



(55) 



with a potential function, V[y%(l), j/2 (0 , • •• ,Un(l)]- The 
RG flow is to the valleys of the potential in the begin- 
ning and then along the valley. The potential structure 
is consistent with the arguments of Ref. 1331 which sug- 
gests that the presence of a potential function is related 
to Zamolodchikov's c-theoremP^ For the commensurate 
case, kj\ — kg — tt/2, the potential is 



V[yi\ = yiaVicV\\c 
- yidy3byu 



■ y\ a y\bVia + yiays, a y\\b - y\ c y\dyi c + yi c y3 C y\\b - yidV\\ a yib 
yicyzbyzc + y\\ c y3 a y3c - yuy3dy\\b + yi a y3 a y3d - y\\ a V3by3d 



,1/2, 2, 2. 2. 2\ ,1/ 2 2 2,2,2,2 2 , 2 \ 

+ 2 {y3a + 2/36 + V3c + V3d + V\\b) VcO + % {^Ic - V2c ~ V\\c + V\\a + Vlb + Via - Vsc + V3d ) Vcir 



(56) 



i„2, 2, 2, 2, 2 I 1^/2 
+ 2 [Vld + Via + Vic + Vlb + y\\b) VsQ + 2 {Vld 



Via 



9 

■Vila 



Via 



Via + Vlb) Vsn, 



where for each Luttinger 

parameter. For Ua 7^ ks or doped cases, we should re- 
move some coupling constants which are not allowed by 
momentum conservation. The RG equations valid even 
when velocities are different are given in Appendix [5] Wc 
checked that both RG equations give consistent results 
when va = vb- 

To get phase diagrams, we integrated the RG equations 
numerically until one of the coupling constant becomes 
of the order of 1(= t). We used initial values of coupling 
constants as small as 10 -8 -10 -3 . Due to the hidden po- 
tential structure, the asymptotic be havi or of the RG flow 
is captured by the following ansatzpEH 

9oi 



9i 



I - L 



(57) 



where l c is the length at which the relevant couplings 
diverge, and the goi determines the ratio among them. 
This represents the fixed ray of relevant coupling con- 
stants. Then, the bosonic fields are pinned down to the 



minima of the effective potential. These values enable us 
to determine the order parameter which takes a non-zero 
value. 

The obtained phase diagrams are shown in Figs. [4]-[8] 
They correspond to the various cases in Table [Tj We also 
investigated the effect of velocity anisotropy to the phases 
in physically relevant parameter region, < J < U/2, 
using the RG equations based on the fermionic Hamilto- 
nian. We studied the range, 1 < va/vb < 10. 

At various points in the following discussion, we also 
label the phases in the "CnS'm" notation, which indicates 
n massless charge modes and m massless spin modes in- 
troduced in Ref. [71 



A. kA — fcs = 7r/2 at half filling 

First, we look at case (a), where fe^ = ks at half fill- 
ing (Fig. El). Absence of the HO phase indicates that 
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CDW 



S-Mott 

-1 

(a)v A /v B = 1 




C1S2 



(b)v A /v B = 5 



FIG. 4. Phase diagram for case (a), k A — fcs = tt/2 with 
equal velocities. Physically relevant region is < J < (7/2. 
A dashed line indicates SU(2) or first order transition. A 
solid line is a Gaussian theory, and a dotted line is an Ising 
transition, "c", "o", "s", and "J" indicate the critical fields 
on the transition line. All of these phases are COST). 



the system does not flow to the enlarged orbital SU(2) 
symmetric state, although this could happen in princi- 
ple by changing initial conditions. The RT phase, a high 
spin sate, resides in < J < U region, which may be 
accessible in real material. 

Although precise boundaries do not coincide exactly, 
the corresponding Hartree-Fock (HF) phase diagram 
shows similar structure.^ There we have the SDW phase 
instead of the RT state in < J < U; they are both 
locally high-spin configurations with anti-fcrromagnetic 
orders along the chain. The S-Mott state found in the 
bosonization result corresponds to mainly the degener- 
ate state of the s'CDW and pCDW orders with a smaller 
region of sSS in the HF phase diagram. In the bosonic 
language, the order parameter of the sSS phase is the 
same as that of the S-Mott state except the total charge 
mode. The s'CDW and pCDW states are degenerate due 
to the orbital symmetry, and they are Ising dual to the S- 
Mott and Haldane orbital (HO) phases respectively (see 
Table IV I . The HO phase does not appear in the bosonic 



calculation, but it is connected to the S-Mott phase by 
the duality mapping of the orbital sector (see Table W\ . 
The CDW state stays almost the same regime in both 
phase diagrams. The BDW phase does not appear in the 
HF phase diagrams, since its HF energy is higher than 
that of the CDW phase. 

The critical properties of the CDW-BDW transition 
are given by a U(l) Gaussian theory of the charge sec- 
tor, while those of the BDW-S-Mott transition are given 
by the Ising theory. The rest of transitions are ei- 
ther SU(2) criticality or first order; the critical fields at 
the BDW-RT transition line are spin modes, and at the 
CDW-S-Mott line are the charge and Ising fields. 

The phase transition from the RT state to the CDW 
state with increasing J is similar to the SDW-CDW tran- 
sition found in the extended Hubbard model (EHM)P^SI 
The EHM has a nearest neighbor interaction, Vrijiij+i, 



in addition to the Hubbard interaction, Unj^nj^. As the 
former interaction becomes predominant, particles try to 
form a CDW state, while strong U prefers SDW. In the 
weak coupling regime, it is found that the SDW state un- 
dergoes a spin-gap transition to a BDW state, and then 
becomes the CDW state through a Gaussian transition of 
the charge sector. In the strong-coupling regime, these 
two transition lines are coupled to a first order transi- 
tion line. In our model, strong J plays the same role as 
V in the EHM; large J induces an attractive on-site in- 
teraction [see Eq. ([2])] leading to the CDW state. The 
properties of transitions from the RT phase to the CDW 
phase, and the existence of the narrow BDW region are 
also the same as in the EHM. Therefore, we expect that 
in the strong-coupling regime, the RT-CDW transition 
in our model also becomes first order, although this has 
not been demonstrated. 

Now, we consider the effect of velo city difference in the 
first quadrant; U, J > [Fig. |4(b)| . As va/vb becomes 
as large as 1.5, we found that the RT and BDW states 
in small J > are completely replaced by a CIS'2 state, 
where only a charge mode of a single band becomes mas- 
sive, and the rest is massless. The CDW and BDW states 
in J > U > are robust to the change in velocity. This is 
because the large anisotropic velocities suppress the in- 
terband scattering, resulting in the domination of intra- 
band scattering. As va/vb is increased beyond 1.5, the 
(715*2 phase becomes larger, although the BDW phase 
always exists between the CDW and CIS'2 phases. 



B. Ica = fcs at incommensurate filling 



CDW 



CDW 



d'SS 



.vSS 



—0 i u °o 

{&)v A /VB = 1 




(b)v A /v B = 5 



FIG. 5. Phase diagram for case (b), k A = ks away from 
half filling with equal velocities. Physically relevant region is 
< J < (7/2. A dashed line indicates SU(2) or first order 
transition. A solid line is a U(l) Gaussian theory, and a dotted 
line is an Ising transition, "c" , "o" , "s" , and "7" indicate the 
critical fields on the transition line. All of these phases are 
CIST). 

The phase diagram for fc^ = ks with incommensurate 
filling is given in Fig. [5] This is similar to the one at half 
filling if we replace the insulating states to correspond- 
ing metallic ones; S-Mott to sSS, RT to d'SS, and BDW 



12 



to CDW. The BDW sate near J ~ in Fig. g] disap- 
pears. The transition between the CDW and d'SS states 
is governed by spin modes leading to the SU(2) criticality 
or first order transition. The CDW-sSS transition is an 
Ising transition, c = 1/2. Finally the d'SS-sSS transition 
is described by a Gaussian theory of the orbital sector. 

The velocity difference in a quadrant, U,J > 0, does 
not modify the large J regime, though a C2S1 state ap- 
pears at small J [Fig. |5(b)| . The C2S1 phase was ob- 
served in other two-leg ladder systems when the velocity 
difference becomes large 

The HF phase diagram of this case 2 ® has the SDW 
phase for < J < Q.6U, which corresponds to the d'SS 
phase found by bosonization; both of them are locally 
high-spin states. In the negative J region, we have the 
sSS state in Fig. [5] while the HF calculation gives not 
only the sSS state, but also a large region of the s'CDW 
phase. As we mentioned, this CDW state is Ising dual to 
the S-Mott phase, which is the insulating analogue of the 
sSS state. In the large J > region, we found p' y -wave 
spin-triplet superconductivity in the HF phase diagram, 
which is replaced by CDW in the bosonic calculation. 
The CDW state around < J < — U is robust, and we 
observe it both at HF level and after renormalization. 



C. h,A 7^ ks at half filling 




S-Mott 



-1 



(! 



1 



FIG. 6. Phase diagram for case (c), kA ks at half filling 
with equal velocities. Physically relevant region is < J < 
17/2. A dashed line indicates SU(2) or first order transition. 
A solid line is a U(l) Gaussian theory. A thin line indicates 
KT transition, "c" , "o" , "s" , and "7" indicate the critical 
fields on the transition line. 

For the system at half filling, but with two differ- 
ent Fermi momenta, the phase diagram is given in Fig. 
|j There is a narrow Luttinger Liquid (LL) phase near 
U ~ J > 0, where all the modes are massless. The transi- 
tion between massive phases and a LL phase is Kosterlitz- 
Thouless (KT) in the sense that a LL phase is critical and 
has power-low correlation, while massive phases have ex- 
ponentially decaying correlations. The transitions among 
Mott phases all have a Gaussian criticality; the S-Mott- 
S'-Mott transition is governed by the charge sector, and 
others are by the orbital sector. On the other hand, the 
s'BDW-D-Mott transition line is the SU(2) criticality or 



first order. In this case, the velocity difference does not 
modify the phase digram in physically relevant region 
essentially. 

The corresponding HF phase diagram shows the SDW 
and s'SDW states in the physically relevant region,®' 
while they are replaced by the s'BDW and D-Mott phases 
in the RG phase diagram. This is a notable difference 
between the result of two equivalent bands and that of 
inequivalent bands. For the Ua = fcs cases, locally high- 
spin states, RT and d'SS, are dominant in < J < U/2, 
while low-spin configurations, s'BDW and D-Mott, are 
found in the kA ^ fc# case. We understand these low- 
spin states as a result of decoherence by two different 
wave numbers. In essence, density waves with different 
phases in the two bands mean that the energy contribu- 
tion from the J interaction averages out to zero. The 
CDW phase, which dominates large positive J region at 
mean-field level, is replaced by the D'-Mott state after 
renormalization. For this case, the pSS state, the metal- 
lic analogue of the D'-Mott state, is subdominant with 
large positive J at the HF level, and this is more en- 
hanced than the CDW order during the renormalization 
flow. Negative J region of the HF phase diagram is again 
covered by the sSS and s'CDW phases, which are related 
to the S(')-Mott state. 



D. kA 7^ fcs at incommensurate filling 



pSS 



.vSS 



s'CDW 



FIG. 7. Phase diagram for case (d), kA 7^ ks away half filling 
with equal velocities. Physically relevant region is < J < 
U/2. A solid line is a U(l) Gaussian theory, and a dotted line 
is an Ising transition. A thin line is KT transition, "c", "o", 
"s", and "J" indicate the critical fields on the transition line. 



With inequivalent Fermi momenta and incommensu- 
rate filling (Fig. [7|, the physically relevant region is cov- 
ered by a LL phase, and the s'CDW phase. The pSS 
and sSS phases can be understood as the reminiscent of 
the D'-Mott and S(')-Mott states which exist at half fill- 
ing. Again, transitions between the LL phase and mas- 
sive phases are KT type except the total charge mode 
remaining massless in both phases. The phase transition 
between the pSS and sSS states is governed by a U(l) 
Gaussian criticality of the orbital sector. The transition 
between the s'CDW and sSS phases is Ising type. As the 
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velocity anisotropy becomes larger, the s'CDW phase is 
gradually suppressed, and whole area in physically rele- 
vant region is covered by the LL state for va/vb > 6. 

The HF phase diagram in this case is similar to the 
RG phase diagram. We have the p(s)SS state with large 
negative U and small positive (negative) J. For large 
negative J, we have the s'CDW state. The pTS and 
SDW states appearing in positive J regime of the HF 
phase diagram are renormalized to the Luttinger liquid 
phase. 



E. kA = it/2 fcs at incommensurate filling 



/;SS 



C1S2 



.vSS 



-1 







FIG. 8. Phase diagram for case (e), kA = 7i"/2 7^ ks away from 
half filling with equal velocities. Physically relevant region 
is < J < U/2. A solid line is a U(l) Gaussian theory, 
and a dotted line is an Ising transition. A thin line is KT 
transition, "c", "o", "s", and "J" indicate the critical fields 
on the transition line. 

The phase diagram of an orbital selective Mott case 
(Fig. [8 1 is almost similar to that of case (d) except Lut- 
tinger iquid state is replaced by a C1S2 phase, where 
a commensurate band opens a charge gap, and the rest 
of the modes remain massless. The velocity anisotropy 
in both directions, va/vb < 1 and va/vb > 1, does not 
modify the C1S2 state in the region, < J < U/2. The 
HF phase diagram for this case is almost the same as the 
previous case with kA 7^ ks at incommensurate filling. 



VII. FULLY SPIN-POLARIZED CASE: 
CHARGE-ORBITAL MODEL 

So far, we have limited our focus to the weak cou- 
pling limit. However, in multiorbital systems the strong- 
coupling limit can bring qualitatively new effectsP^H^ 
In particular, in the physical (J > 0) case, one may ex- 
pect every ion to be in the state of maximal spin consis- 
tent with given total occupation. In this circumstance, 
we expect antiferromagnetic order in the half-filled case, 
and when the system is slightly doped, more complicated 
structures such as phase separation and spiral phases will 
appear. At filling further away from n = 2, the system 
shows ferromagnetism (FM) with an orbital order. Con- 
sidering the fact that ferromagnetic states dominate the 



phase diagram at most fillings ESEJU i n this section, we 
investigate the possible orbital orders assuming that the 
system is fully spin polarized. In other words, we con- 
sider the effect of residual backscattcring in the subspace 
of the charge and orbital sectors assuming the spin ex- 
citations are frozen. We leave the investigation of the 
regime close to half filling for future study. 

Suppose that all the electrons have the same spins, the 
model in Eq. ([I]) and Eq. |2]) is then reduced to, 



H = ^kinetic + z2(U- 3 J) n iA n iBl 



(58) 



where we omit the spin index. Now the SU(2) S symme- 
try is lost, but we can regard the orbital part as pseudo 
spins. If the two bands have the same kp and velocities, 
the system has the orbital SU(2) symmetry. The band 
splitting in orbital sector is isomorphic to the Zeeman 
splitting by magnetic field. Therefore, the model now 
turns to a simple Hubbard model with effective interac- 
tion U e s — U — 3 J, with or without magnetic field. 

Depending on the effective crystal field splitting be- 
tween the two orbitals A, and the band widths, there 
may be three different scenarios in this model. The first 
case is that the two bands are completely degenerate: 
A = 0, and fc^ = ks- We expect a staggered orbital 
order to appear. On the other hand, when there exists 
either small splitting or when the band widths are slightly 
different, the two momenta are not equal, and orbital or- 
ders might be suppressed. We will discuss these scenarios 
below using bosonization method. However, there is an- 
other scenario, which may arise when either splitting is 
large, or two band widths are greatly different. Then, 
only one band has states at the Fermi surface and the 
physics is trivial. 

We first look at the degenerate case, kA = fcs, corre- 
sponding to the absence of magnetic field. The bosonized 
form of the Hamiltonian in Eq. ( 58 1 at half filling is given 
by 



H 



K v {ve v f + -L(v<^) 2 



+ (2^ C/cffCOS ( 2 ^°) 
■ U c s cos (2%/20 c ) 



(59) 



(27ra) 2 

where the Luttinger parameters and velocities are, 



v c (o)/K, 



K c v c = K v = v 



v\l± 



(60) 
(61) 



Thus, the charge and orbital modes are decoupled, and 
each mode has a SU(2) symmetry. The total symmetry 
is SU(2) x SU(2) = SO(4). 

Translating the analysis for the Hubbard modeP2 to 
our charge-orbital model, we found the following results. 
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orbital sector is always massless. With an attractive in- 
teraction, the band degeneracy occurs with smaller fill- 
ing than with a repulsive interaction, since we guess the 
upper band is pulled down by the lower filled band for 
U c s < 0. The charge mode is massive (massless) for 
repulsive (attractive) interaction at half filling. At tree- 
level, the states with longest correlations are ODW x and 
ODW y for U cS > 0, and p' x TS and p' y TS for U eS < 0. 
Since the orbital symmetry is explicitly broken, the ex- 
ponents of correlation may differ for different directions. 
Finally, contrary to the complete degenerate case, we 
speculate that an orbital selective Mott phase appears 
near n = 1 for inequivalent bands; once one band is half 
filled, the commensurate wave vector opens a gap, and 
the other band remains metallic. Further filling just goes 
to the metallic band until it reaches half filling. Fig. |9(b)| 
presents the general phase diagram for this case. 



FIG. 9. Schematic phase diagrams for the spin-polarized 
model for two equivalent bands or slightly different bands. 
When h,A 7^ ks , the ODW^ and ODW H phases have the same 
correlation exponent, and so do p' x TS and p^TS. Filling closer 
to n — 2 is not investigated. "CnOm " represents a state with 
n massless charge modes and m massless orbital modes. 



For incommensurate filling, the last Umklapp term in Eq. 
(591 vanishes, and the charge mode is massless. Also the 
SU(2) C symmetry is broken to U(l) c . About the orbital 



sector, we find: 

(1) Ues > 0: Orbital density wave (ODW) has the 

longest correlation, and both orbital and charge 
modes are massless. The SU(2) D symmetry requires 
ODW about all three directions (x, y, z) are degen- 
erate. 

(2) f/ e ff < 0: The orbital sector becomes massive, and 

the phase with slowest decaying correlation is 
orbital-singlet superconductivity with parallel spins 
i.e. p' y TS in Ref. EH 

At half filling, the charge mode becomes massive (K c = 
1/2) when the effective interaction is repulsive; the sys- 
tem is Mott insulating. The orbital part still gives ODW, 
and this FM+OD W sta te in U > 3 J re gime is observed 
both analyticall y 28 1 47 1 481 and numerically. 1 ^ 2 ' For the at- 
tractive side, the charge mode is gapless [K c = 1) with 
an orbital gap by the P^TS order; this is the Luther- 
Emery phase. This triplet superconductivity agrees with 
the numerical result by Sakamoto et aZ.p3 The results 
are summarized in Fig. 



9(a) 



Now we turn to the case with Ua ^ kg; there is a 
pseudo magnetic field acting on the orbital space. At 
very small filling, only a single band is filled, so the 
ground state is ferromagnetic Luttinger liquid of a sin- 
gle gapless mode. When we dope enough the two bands 
start to share the Fermi surface. The SU(2) Q symme- 
try is reduced to U(l) , and the cos (2-\/20 o ) term van- 
ishes due to two different Fermi momenta. Thus, the 



VIII. CONCLUSION 

In this paper, we investigated a two-orbital Hubbard 
model which may encapsulate phenomena realized in 
transition-metal nanowires. Along with many aspects 
of two-leg ladder models, unique properties of transition- 
metal d-orbitals lead to several new results. 

In our analys is, we us ed the ideas of dynamical symme- 
try enlargement! 10 ' 20 '^ and duality relation ^ 22 ' 23 ! to list 8 
insulating phases at half filling, and 6 metallic phases 
away from half filling. Each phase is represented by a 
Gross-Neveu model at low energy, and phases are related 
to each other by duality mappings. The same analysis 
was done for weak- tunne ling models with the U(1) D sym- 
metry about z- axis J2H2U anc | f or strong-tunneling models 

with the U(1) Q symmetry! 10 ' 16 ' 17 '. Some cases we stud- 
ied in this paper belong to the latter category, and we 
obtained the same results. However, we also studied a 
model with a different orbital symmetry, the U(1) Q sym- 
metry about y-axis, and found that this form of inter- 
action gives rise to a new combination of ground states; 
half of which were found in weak-tunneling systems and 
the other half were in strong tunneling systems. This is 
because the underlying orbital rotational symmetry in- 
volves a generator which cannot be written by a single lo- 
cal bosonic field. The quantum phase transitions among 
those ground states were also briefly summarized using 
Gross-Neveu models as low-energy effective theories! 16 ' 22 ! 
The new combinations of ground states we found means 
that the model we proposed can exhibit new kinds of 
quantum phase transitions. 

After classifying the possible ground states, we deter- 
mined the ground state of the model by numerically inte- 
grating the RG equations. The phase diagrams are pre- 
sented in Figs. |3^[8] The essential results for physically 
relevant regime, < J < U/2, are as follows: 

(1): For two equivalent bands at half filling, we found a 
high spin state, RT, for relatively large J, and a 
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low spin state, BDW, for small J (Fig. [IJ. Keep- 
ing k,A = ks but introducing a velocity difference 
causes these phases to be replaced by a C1S2 state 
where only a charge mode in a single band becomes 
massive. Upon doping, the RT state turns into 
the d'SS state and the BDW state disappears (Fig. 
[5]). A velocity difference leads to a C2S1 phase for 
small J, but the d'SS state still survives in large J. 

(2): For the cases with two different Fer mi mome nta, 
similar RG analyses have been doneP^ 16 * 17 ! Our 
main contribution for this case is that we identify 
the ground states for our specific model and in- 
teractions. The ground state in the physically rel- 
evant region is the s'BDW state for large J and 
the D-Mott state for small J at half filling (Fig. 
[6|. These phases are robust to velocity differences. 
The surprising result here is that we have locally 
low-spin configurations even when J > 0, while 
completely degenerate bands give high-spin states, 
RT and d'SS, in the same parameter regime. We 
think the low-spin configurations are achieved by 
the destructive interference between two different 
Fermi momenta, and therefore, the complete or- 
bital degeneracy is crucial to have locally high-spin 
states. When the system is away from half fill- 
ing, a Luttingcr liquid phase becomes dominant, 
although we observed the s'CDW phase for very 
small J (Fig. FFb. This density wave state is wiped 
away by velocity difference, and replaced by the 
Luttinger liquid state. When only a single band 
is commensurate, the system exhibits an orbital- 
selective charge-gapped state due to the intraband 
Umklapp process while all the other modes remain 
massless (Fig. [8]). 

We also investigated the charge-orbital model obtained 
as a result of full spin-polarization. The model can be 
mapped to a Hubbard model with effective interaction 
U eS = U- 3 J. For repulsive side, we find the FM+ODW 
phase with massless charge and orbital modes is stable 
at most of the fillings. When two bands are equivalent, 
we have a charge gap at n = 1, while we expect an or- 
bital selective Mott phase around n = 1 for inequiva- 
lent bands. The attractive side is dominated by orbital 
singlet superconductivity with an orbital gap when two 
bands are equal. On the other hand, for inequivalent 
bands, orbital singlet and triplet superconductivity have 
the longest correlation and the orbital gap is absent. The 
charge gap does not develop in the attractive side. 

Finally we refer to possible experimental consequences. 
From the band calculation, it is plausible that the real 
system exhibits a ferromagnetic state, meaning that the 
system has relatively strong correlations!^ In particu- 
lar, since U is much larger than J in real materials, we 
guess the best possible ground state is a ferromagnetic or- 
bital density wave from the discussion in Section |VII| Of 
course, we have to note that the result is based on strong 
assumptions such that the system is fully spin polarized, 



and that the number of the bands is just two. 
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Appendix A: Abelian Bosonization 

In this appendix, we briefly recall formulas from 
Abelian-bosonization which are needed in the textP^U 
The following formula gives transformation from a 
fermionic Hamiltonian to a bosonic Hamiltonian 

i, mar = Jj^e 1 **— (Al) 
V 2ira 

where m — A, B is orbital, and r = i?, L is chirality. The 
bosonic fields satisfy the commutation relations, 

[$m<rR(L)(x), ^ m 'o'R(L){x')] = ±Mr<S mm '<Wsgn(a; - x') 

(A2) 

The Majorana fermions take care of the fermionic prop- 
erties and obey anti-commutation relations 

\y]m<7 : Tim' a' } — 2S mrn ' &ao' • (-^^0 

A more convenient representation is given by the non- 
chiral fields, 

4>ma,0 ma = ^(<I> m(T L ± $ m<7 fl). (A4) 

They are connected to density and current as V</> oc n 
and V# oc j, and satisfy commutation relations, 

[4>ma(x),0 m 'a'(x')] = -iTtS mm >S a(T >Q(x' - x). 

Finally we move to different combinations of these fields, 



4>c0 




"l 1 1 1 " 




<t>A\ 


4>C7T 


1 


11-1-1 




<t>M 


4>s0 


_ 2 


1-11-1 




<f>Bt 






1-1-1 1 







(A6) 



where /i = (c, s) represents charge or spin modes, and 
v = (0, 7r) give bonding/antibonding basis. #'s are trans- 
formed in the same manner. 

The sign of each coupling constant is determined by 
Klein factors, and by a commutator between different 
chirality, [<f>R >n (x), <&L,n'(x')] = i^^n.n'- The eigenvalues 
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of Klein factors composed of two Majorana fermions (dif- 
ferent from the ones introduced for refermionization) are 
taken to be 



% = VAsVBs = VAtVAi 

= VAtVBl = VBtVAl = -VBtVBi- (A7) 



Appendix B: RG equations for unequal velocities 



When two Fermi velocities are different, it is more con- 
venient to use current operators than using refermion- 
ization. We follow the notation of Ref. 7 with slight 
modification: 



Jrar — ^ ^msr^ms' r 5 J mr — ^ ^rnsr® ss'lpms'r (Bl) 

ss' ss' 

L r = /,j> Asr i>B s >r, L r = ^ j)\ sr (T ss 'lj)Bs'r (B2) 

ss' ss' 

M mr = -ilprntr^mlr, N rss > = ^ rAs 1prBs' ■ (B3) 



When fc/i ^ ks, the interactions terms are given by, 

— %int = QlpJ ArJ AL + giaJ AR ' J AL 
+ 92 p JbrJbL + 92a J BR' J BL 

+ 9x P {JarJbl + JbrJal) (B4) 

+ 9xa {J AR ■ J BL + J BR ■ J Al) 

+ g t p(L R L L + h.c.) + g t a(L R ■ L L + h.c). 

This expression is formally the same as the one given in 
Ref. GO When k A = fcs, we have additional processes, 



g ap (l r l\ + h.c.) + g aa ( 



Lr ■ L\ + h.C, 



(B5) 

Umklapp processes are allowed when the filling is com- 
mensurate (n = 2): 



-Kt = 9m (MarM al + h.c) 
+ 92u (m B rMI l + h.c 



[M ar mI l + M BR M\ L + h.c. 



(B6) 



+ 9tu P (N Rap N Lafs - N Ra/3 N L p a + h.c 

+ 9tua (rfiapNLaf) + N^pN^a + h.C.) . 

The ~g\ u and ~gi u processes are allowed only when each 
band has commensurate filling, i.e., k m = tt/2. We ig- 
nore all the chiral scattering processes, since they only 
renormalize the velocities. 

In the following, we use the renormalized coupling con- 
stants, yi = g~iiT~ 1 (vA+VB)~ 1 ■ The RG equations for the 
k A = ks = 7r/2 case are 



yip = -P {vlp + 3yL + ^vLa + yLp - yip - 3 2/L) - a viu 

V2 P = -a (y 2 + 3yl a + iy 2 tua + y\ 



Vxp — y ap + 3y QCT 3y tua y tup y tp 



Vtp 
- 3yl 



3yia) - Pvlu 



Vxu 



(B7) 



y\cT = -2/3 (y aa (y aa + y ap ) + y tua (y tucr + y tup ) + y ta [y ta - y tp )) - 4ay 1(T 

ho- = 

Vxa = 



-2a (y aa (y aa + y ap ) + ytua (ytua + ytup) + yta (yta - yt P )) - 4ft/f 

"2 (Vaa (Vaa ~ Vap) + Vtua (ytua ~ ytup) + Vta (yta + Vtp)) ~ 4y, 



2 

xa 



(B8) 



ytp = -2ytu P yxu + yt P y c - + iyt a y s - 

yta = 2y tua y xu + y tp y s - + yta (y c - - 2y s+ ) 

y ap = -ytup (ayiu + Py2u) - y ap y c - - 3y aa y s - 

llaa = -ytua (otyiu + PV2u) - VapVs- - Vaa (Vc- + tys+) 



(B9) 
(BIO) 



yiu = -4 (3(3y a aytua + Py apytup ~t~ Qyiuyip) 
i)2u = -4 (Sayaaytua + ay ap ytu P + fiy2 U y2p) 

9xu — 4 (3y ta ytua — ytpVtup — yxuVxp) 



(Bll) 
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Vtup = -y ap {ayiu + Pv2u) - ZytpUxu - ytu P y c + - 3y tticr j/ s - 

ytua = -y a <r (ayiu + fiy2u) + ty tcr y xu - y tup ys- - ytua (y c + + 2y s+ ) 



(B12) 



where we defined y c ( s )± = ay lp{<j) +f3y 2p (a) ±2y xp (a) with 
a = (v A + Vb)/(2va), and fi = (v A + v B )/(2v B )- For 
doped cases, and Ua ^ k B cases, the coupling constants 
which are not allowed by momentum conservation should 
be removed. 

As we mentioned, the asymptotic behavior of a RG 
flow is captured by the ansatz (57), and now the ratios 



of coupling constants at fixed points depend on velocity 



r 



differences. However, we can easily distinguish phases 
with different fixed point structure by looking at the signs 
of relevant couplings, and irrelevant couplings. In that 
sense, we identify phases as the same ones when the rel- 
evant couplings and the signs are the same. When the 
relevant couplings are different, or the signs of renormal- 
ized couplings are different, we regard them as different 
phases. 
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